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Abstract 

In this work, we review and clarify computational issues about the VK-gaugc boson one- 
loop contribution to the if — )■ 77 decay amplitude, in the unitary gauge and in the Standard 
Model. We find that highly divergent integrals depend upon the choice of shifting momenta 
with arbitrary vectors. One particular combination of these arbitrary vectors reduces the 
superficial divergency down to a logarithmic one. The remaining ambiguity is then fixed by 
exploiting gauge invariance and the Goldstone Boson Equivalence Theorem. Our method is 
strictly realised in four-dimensions. The result for the amplitude agrees with the "famous" 
one obtained using dimensional regularisation (DR) in the limit d — > 4, where d is the number 
of spatial dimensions in Euclidean space. At the exact equality d — 4, a, three-sphere surface 
term appears that renders the Ward Identities and the equivalence theorem inconsistent. 
We also examined a recently proposed four-dimensional regularisation scheme and found 
agreement with the DR outcome. 
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Figure 1: VF-gauge boson contribution to the H 
with relevant shift vectors are indicated. 



77 amphtude. Momentum flow together 



1 Introduction 



Today one of the main focal points at the Large Hadron Collider (LHC) is to search for the 
Higgs boson (H) fl[p' 
Indeed, the recent J&Jj 



through its decay into two photons, — )• 77 (for reviews see [4|[5]). 
observation by ATLAS and CMS experiments of a resonance, that could 
be the Standard Model (SM) Higgs boson, is based on data mainly driven by — >• 77. In the 
(SM) (8 -10 , this particular decay process goes through loop induced diagrams involving either 
charged fermions or M^^-gauge bosons. Their calculation was first performed in ref. 11 in the 



limit of light Higgs mass mu <C using dimensional regularisation in the 't Hooft-Feynman 
gauge. Since then, there are numerous works spent on improving this calculation including finite 



Higgs mass effects in linear and non-linear gauges 12 -14 , different regularisation schemes 15 -18 
and/or different gauge choices flo]. 

The -?/—)• 77 amplitude is originated, in broken (unbroken) phase, by a dimension-5 (dimension- 
6) SM gauge invariant operator(s) and, therefore, its expression, within a renormalizable theory, 
must be finite, gauge invariant and independent of any gauge choice. The amplitude should also 
be consistent with the Goldstone Boson Equivalence Theorem (GBET) [20| - [22| since the SM is 
a renormalizable, spontaneously broken, gauge field theory. 

A problem arises when the Pi^-gauge boson contribution (see Fig. [T]) to -ff — )• 77 produces 
"infinite" results at intermediate steps. These problems are usually treated by using a gauge 
invariant regulator method, e.g., dimensional regularization. In the unitary gauge 23 , this 



indeterminacy is more pronounced and more difficult to handle with due to the particular form 
of the VF-gauge boson propagator. On the other hand it is much simpler to work with only few 
diagrams, that involve physical particle masses, rather than many. 

More specifically, in the unitary gauge, one encounters divergencies up to the sixth power. 
It is well known that, in four-dimensions, shifting momenta in integrals that are more than log- 



2 



arithmically divergent is a "tricky business" - recall the calculation of linearly divergent fermion 



triangles in chiral anomalies 24 , 25 - that requires keeping track of several "surface" terms for 
these integrals. There is also the situation we face here where apparent logarithmically divergent 
integrals turn out to be finite but discontinuous at d = 4. 

We would like to bypass those ambiguities and at the same time to present a "regularisation" 
method, by performing the calculation for the — )• 77 amplitude strictly in 4-dimensions and 
in the physical unitary gauge. Our method is similar to the one used elsewhere for calculating 
triple gauge boson amplitudes [26j[27j , or Lorentz non-invariant amplitudes [28j , and consists of 
three steps: 



1. We write down the most general Lorentz invariant — )• 77 amplitude. 

2. We introduce arbitrary vectors that account for the "shifting momentum" indeterminacy. 
We show that a particular choice of those "shifting vectors" cancel higher powers of in- 
finities leaving still behind at most logarithmically divergent integrals that are treated as 
undetermined variables. 

3. We exploit physics, i.e., gauge invariance (Ward Identities) and the GBET in order to fix 
the last undetermined variables. 



This method is quite general and can be applied to other observables too. Following these steps 
we arrive at the same result for the — t- 77 amplitude obtained by J. Ellis et.al and by 
M. Shifman et.al [l3] almost 35 years ago. Our analysis, among other issues, highlights that 
the recent observation (6|[7] of the iJ — )• 77 at the LHC signifies the validity of the Goldstone 
Boson Equivalence Theorem. As a further clarification we also make a remark on the direct 
calculation in the following three cases: we first perform the integrals in exactly d = 4 (with 
no regularisation method beyond the one discussed in point 2 above), second, by exploiting 



Dimensional Regularisation (DR) as defined in refs. 29 ,30 and then taking the limit d — >• 4, and 



finally third by using a four-dimensional regularization scheme introduced in ref. 131 



Our calculation is complementary to, but somewhat different than, the two existing ones 19 



32 , 33 performed in the unitary gauge. Our results agree with ref. 19 . 

The outline of the paper is as following: in section [2] we present the calculation of the VF-loop 
contributioi:[^ to — )• 77 amplitude, its ambiguities and the resolution within physics arising 
from GBET. Next, in section [3] we examine details of the amplitude calculation within an alter- 
native and recently proposed four dimensional regularization scheme [31] , the one that resembles 
most closely the symmetry approach taken here. In section |4] we discuss other possible physical 
setups plus experiment that may help to resolve inconsistencies. We conclude in section [5j There 
are two Appendices : A that contains some intermediate formulae and B where we present the 
details about the origin of surface terms in four dimensions. 



^Note that the calculation of the fermion triangle contribution is well defined i.e., it is independent of arbitrary 
vectors and finite. We are not going to repeat this calculation here and refer the reader to the reviews cited above. 
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2 The VF-loop contribution to i7 — >^ 77 in SM 

The most general, Lorentz and CP- invariant, form of the of-shell — )• 77 amphtude is, 

Mig^" + M2 K ^2 + -^3 < + Ma A:f k'( + M5 ki^ k!^ , (2.1) 

where ki and k2 are the outgoing photon momenta shown in Fig.[T| and the coefficients A^j=i..5 = 
Aii=i,,5{ki, k2) are scalar functions of k1,k2, and ki ■ k2. By considering that all particles are 
on-mass-shell, that is A;^ = = , fei • ^2 = '71^/2 , ki ■ e*{ki) = ,k2 ■ e*(fe2) = 0, we obtain an 
amplitude Ai = A^'^'^e* (A:i)e* (^2) with only two, undetermined (for the time being), coefficients, 

M^"' = Mig^"' + M2k'(k^ . (2.2) 

In unitary gauge, the Feynman diagrams that contribute to Mi and M2 are displayed in Fig.[T] 
In order to calculate them, we introduce three arbitrary four- vectors a, b and c, one for each 
diagram. These vectors shift the integration momentum, i.e., p ^ p + a for the ffist diagram, 
p p + b for the second diagram and p ^ p + c for the third diagram. As we shall see, these 
arbitrary vectors operate as regulators capable to handle highly divergent integrals related to 
unitary gauge choice. Furthermore, the vectors a, b and c, linearly depend upon the external 



momenta ki and ^2. Hence a,b and c are not linearly independent [c.f. eq. (2.4)]. This is an 
important fact leading to the cancellation of infinities. 



We first calculate the less divergent part of M'^'^ in eq. (2.2) which is the M2 coefficienlj^ 



By naive power counting, we see that M2 diverges by at most four powers. Then we perform 
the Feynman integral calculations strictly in 4-dimensions. For reasons that will become clear 
later, we shall keep the number of dimensions general in all intermediate steps of the calculation 
i.e., g^^g^u = d. As we will see, d contributes only in finite pieces of M2 [c.f. eq. (|2^g 
With all the above definitions, we can write down the total amplitude in the form 



+ A21 {p + a)'' {p + aY + A22 {p + bY {p + by + ^23 {p + cY {p + cY 

+ ^31 {P + aY k'i + A^2 {p + bY k'( + ^33 {P + cY k'i 

+ Aai {p + aY k!^ + Aa2 {p + bY k!^ + ^43 {p + cY kl^ 

+ Aifc^^r] , (2.3) 



where the coefficients Aij = Aij{p'^;ki,k2;a;b;c) with —6 < n < 0, are given in Appendix A 
and the ~ sign is the proportionality factor: Note that M^'^ is a (superficially) 6th power 



divergent amplitude in the unitary gauge, in eq. (2.3) solely contributes to Mi while all 



other ^-elements contribute to both Mi and/or M2 in eq. (2.2). 



First we focus on the calculation of the "less divergent" coefficient M2 of eq. (2.2). Based on 



naive power counting, we observe that the ^21, -^.22, -4.23-terms in eq. (2.3), lead to at the most 

^The coefficient Mi will be fixed later on by the requirement of gauge invariance. 

^On the contrary, we shall see that there are non-trivial d-contributions into A^i -coefficient. 
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quartic divergent integrals. However, when adding all these pieces together, we find that quartic 
divergent integrals vanish for every arbitrary vectors a, b and c leaving behind an expression 
with integrals of third power (in momenta) plus integrals with smaller divergencies. Then the 
cubically divergent integrals are proportional to all possible Lorentz invariant combinations like: 
[{a + b- 2c) • p] p^p^ , [{a + b- 2cYp'^] p^ and [(a + 6 - 2cYp^] p'^. Therefore, choosing 



a + 6 - 2c 



0, 



(2.4) 



we ensure that third order divergent integrals related to ^21) -^22, -423-terms, vanish identically. 
In the same way, by naive power counting, ^31 and ^33-terms - these terms in eq. (2.3) together 
with A32 contribute solely to 7W2 - lead again to at most third order divergent integrals. How- 



ever, in the sum of ^31 and ^33-terms in eq. (2.3), third order divergent integrals vanish for 



arbitrary a, b and c, leading to an expression, that when added to ^32-term, consists of at most 
quadratically divergent integrals, proportional to [(c — a) ■p]p'^k'( and [(c — a)^/cj']p^. We choose. 



0, 



(2.5) 



for the quadratically divergent integrals to vanish. Likewise, when we add A42 and ^43-terms - 
these terms, together with A41, solely contribute to AI2 in eq. (2.2) - the third order divergent 



integrals vanish for every choice of a,b,c leading to an expression, that when added to A41, 
consists of at most quadratically divergent integrals proportional to [(c — b) ■ plp'^kl^ and [(c — 
bykl^jp'^. Therefore, we choose 

c-b= , (2.6) 



for infinities to vanish identically. Prom eqs. (2.4), (2.5) and (2.6) we arrive at the final relation 
among the three introduced vectors: 

a = b = c. (2.7) 



Eq. (2.7) suggests that the rest of the divergent integrals depend by, at most, one arbitrary 
vector, say the a- vector. Note that ^51 contributes only to the finite part of A^2- Now, if we 



impose conditions (2.7) onto the remaining expressions for A2i,A22, -Asi-terms of eq. (2.3), we 



find that all quadratically and linearly divergent integrals vanish, independently of the direction 
of the a-vector. We stress here the fact that the cancellation of divergencies down to logarithmic 
ones is a highly non-trivial, almost "miraculous", result. These cancellations only take place 
for a particular choice of the momentum-variable shift vectors, [eq. {2. 7\ )]. Of course this is an 
expected outcome for an observable in a renormalizable theory. 

The final result contains at most logarithmically divergent integrals. Despite of the fact that 
the resulting expressions so far contain the shift j? + a instead of p with an arbitrary vector a, 
its presence is irrelevant since logarithmically divergent integrals are momentum-variable shift 
independent 34 . This result is different with the one obtained in refs. |32||33| , where there 
is a quadratically divergent term remaining and is tuned to zero by appropriate choice of loop 
momentum. Summing up all the above contributions to M2, we find a particularly nice and 
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symmetric form for Ai^'^ , 



M 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



(2^) 



4((i- + 2mf^ 



4((i-l) m^ + 2mf^ 



[p2 - mly] [{p - - m^] [{p -ki- k2Y - m 



2 1 



d'^p 



(27r)4 P'^^'^ \ [p2 _ ml,][{p - fci)2 - ml,][ip - - k^? - m'^] 

-4 {p ■ k2) 

[p2 - m2^] [{p - k2Y - m^] [{p -ki- k2Y - m"^] 



-4(d - 1) - 4 (p • A:2) 



d^p ^ 



(2^) 



-A{p-ki) 



p^k^ 

4 2 1 _ - A;i)2 - m^][(p - A;i - A;2)2 - m? 



A{d-l)ml,-A{p-ki] 



[p2 - rn^] [{p - ^2)2 - m^] [{p - ki - /c2)2 - m 



d'^p 
(2^ 



^1^2 



6mL + 2p2 



[p2 - - A;i)2 - m^][(j< - /ci - A;2)2 - rn^] 



6mL + 2p2 



[p2 - 777,2 ] [(-p _ ;j2)2 _ 7,72 ] Jj-^, _ /j^ _ /j2)2 - m' 



(2i 



Introducing Feynman parameters, shifting momentum variable from pto I and ignoring all terms 
that contribute to A^i |^we find that the contribution to ^A2 in eq. (2.2) arises solely from the 
term, 



M2klki^ ~ 8 / dx I dy 
Jo Jo 



dH f k'(k^-2{£- k2) £t'k'(-2{i- ki) r k^ 



(£2 - A)= 



.1 rl-x 
+ 8 / dx 

Jo JO 



dy 



2{d-l)x{l-x-y) 



k'(k>^, (2.9) 



with A = x{x + y — l)rn'jj + Obviously, the first integral in eq. (2.9) is (superficially) 

logarithmically divergent while the second one is finite. The number of dimensions (d) appears 
only at the finite integral and therefore we can fearlessly set d = 4 everywhere. This means that 
we do not use dimensional regularisation in what follows (see however the discussion below). We 
state here few additional remarks to be exploited later on: a) we observe that the top line in 
the integrand of eq. (2.9) does not vanish in the limit — )• and, b) despite of appearances 



in eq. (2.8), there is no mjj in the numerators of the subsequent expression eq. (2.9). The whole 



m'jj contribution arises from the denominator's A-term. 



These terms will be used later in arriving at eq. (2.241 
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Our next step is to parametrize the logarithmically divergent integral in eq. (2.9) by an 
unknown, dimensionless, parameter A to be determined later by a physical argument. So we 
define, 



dx 



dy 



d^e £2 k'iki^-2{l- k2) £^'k'[-2{£■ ki) t k^ _ iX 



4(4^) 



k\ k^ . 



(2.10) 



An important parenthesis here. We could of course promote d'^l — )■ d'^t and use dimensional 

—7- \(^g^^ in d-dimensions. In this 
-1 which is finite and non-zero, and, 



regularisation 



29 



by exploiting symmetric integration 
case, and after taking the limit d — )• 4, one finds A = 



agrees with the one we find below in eq. (|2.20|) after imposing the GBET condition. This is also 

the integral 



However, according to refs. [32 33 



the result found in the original refs. [IT- 13 

in eq. (2.10) is discontinuous at d = 4; in fact, when symmetric integration 
d 



m 



4 is used, one finds instead A 



0. This is also understood in a slightly different context. It 



has long been known 134 - 36 that shifts of integration variables in linearly (and above) divergent 



integrals are accompanied by "surface" terms that appear only in four dimensions - a famous 
example being the integrals in chiral anomaly triangle graphs. For our purpose here lets start 



with the following shift of variables in a linearly divergent integral that has been generalised 36 
to work in 2(jj-dimensions following the expression, 



[(^-A;)2- A]^ 



(£2 - A) 2 



ku (5, 



(2.11) 



that is valid for < 5/2 and A constant, possibly dependent on Feynman parameters, like the 



one given below eq. (2.9), a nd A:^ is an arbitrary constant four vector. By taking the derivative. 



of both sides in eq. (2.11) and shifting the integration variable for the logarithmically 



divergent integral encountered, and evaluating the finite one, we easily arrive at 



ITT 



TT 



uj-2 



r(3 



(2.12) 



For an alternative and detailed proof of eq. (2.12), see Appendix Bp Applying eq. (2.12) to i'^£P 
terms of eq. (2.10) with (^^! ^^d. 



0, 



UJ 



-1, uj = 2-e (DR) 



(2.13) 



This is consistent with the symmetric integration in 4-dimensions (cj = 2), but, is also consistent 
with the usual tabulated textbook result [37] from dimensional regularisation in 4— 2e-dimensions 
Eq. (2.13) shows that A is discontinuous at d = 2tt; = 4. Then the Question arises: 
the one that is indicated by well defined, calculable, boundary 



e . 



which A to believe in? Answer: 
conditions and symmetries of the underlying theory. 

^The same result is obtained by standard algebraic tricks. We would like to thank R. Jackiw for communicating 
his calculation to us. 
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The above parenthesis to our calculation motivates us to avoid the direct calculation of 
integral (2.10) but set d = 4 everywhere and treat A as an unknown parameter to be defined 



later within a physical context or experiment. Substituting eq. (2.10) into eq. (2.9) we arrive at 



87r2 



A — 6 m 



w 



dx 



dy 



1 - 22;(1 



y) 



(2.14) 



Evaluating the double finite integral in eq. (2.14), and restoring the proportionality factor given 



below eq. (2.3), we obtain. 



where 



/3 



4m 



W 



m 



and, /(/3) 



-2A + 



3/3 + 3/3(2-/3)/(/3) 



arctan^ ^-^^=^ , f3 > I 



H 



In 



^ vr 



(2.15) 



(2.16) 



/3 < 1 



Our final step is to determine the unknown parameter A in eq. (2.15). For this we need 
physics that reproduces A^2 in a different and unambiguous way. One choice is to adopt the 



Goldstone Boson Equivalence Theorem (GBET) [20-22 which states that the amplitude for 
emission or absorption of a longitudinally polarised W at high energy becomes equivalent to the 
emission or absorption of the Goldstone boson that was eaten. Mathematically, this is written 
by an equation [38j , 



S[W^, physical] =es[s^, physical] 



(2.17) 



which says that the S-matrix elements for the scattering of the physical longitudinal vector bosons 
Wl with other physical particles are the same as the S-matrix elements of the theory where the 
Wl's have been replaced by physical Goldstone bosons (s^). We are not going to get into details 
apart from the original literature the reader is also referred to the articles [19 

771' 



38-40 



w 



Is 



0, 
0, or 



here 

Following ref. [38j, within perturbation theory and in the limit of high energies 
GBET can be expressed with physics in two different limits of the theory: (a) q^/Xh 
(b) m]j/s 0. 

The limit (b) is irrelevant]^ for defining A in eq. (2.15) so we completely focus on the limit 
(a). It is very easy to see that, in the unitary gauge, the Wis do not decoupl^for vanishing 
gauge coupling g. Consider for example the diagrams in Fig. [ij there is always a from the 
HWW-veitex. that cancels another m'^r sitting in the denominator of the longitudinal part for 
the internal W-boson propagator expression written in the unitary gauge. So, as it was already 
noted in the paragraph below eq. (2.9), in the limit 5 — )• there are remaining non-decoupled 



terms. Unfortunately, these effects may be obscured or misjudged by the regularisation method 
needed to handle divergent, intermediate, loop integrals. This is exactly what happens here 



^The limit (b) simply says that matrix elements for the theory which contains the physical Wl's and zero v.e.v 
is equal to those produced by scattering of massless physical Goldstone bosons (instead of Wl's) at high energies. 
We have checked that eq. (2.171 is satisfied in this limit. 

''This is another advantage of calculating in the unitary gauge. 
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Figure 2: Charged Goldstone boson contributions to i/ — t- 77 in the Hmit of g — >• 0. 



when trying to calculate A directly from its ambiguous form (2.10). On the other hand however, 
at the exact g = 0, with fixed v.e.v v and Higgs quartic coupling A^, eq. (2.17) suggests that 
the Goldstone bosons (s^) should reappear at the physical spectrum of the theory while the 
longitudinal components of W^s become unphysical. At this limit, the SM is a spontaneously 
broken global SU{2)l x C/(l)y-symmetry that couples, minimally, to electromagnetism. The 
interactions between the Higgs and photon with the Goldstone bosons are simply those of a 
spontaneously broken scalar QED with [/(l)om, 



Hs+s- 



im 



H 



7s+(pi)s (P2) : -ie{pi+P2Y 



77 s s : 2ie g 



2 uv 



(2.18) 



Armed with these Feynman rules we calculate the diagrams in Fig. [2] By doing so, we 
introduce again three momentum variable shift vectors, one for each diagram, exactly in the 
same way we did for the calculation of the diagrams in Fig. [T] The Lorentz structure of the 



amplitude is completely analogous to eq. (2.2) with A^i,2 -^i,2(gbet)) but now due to the 



scalar propagators, the superficial degree of divergence, for diagrams contributing to A^2(gbet)) 
is -D = —2. Hence, all integrals involved in A^2(gbet) ^-^^ finite and in addition, they are 
independent of any momentum integration shift vector variable. As a consequence, A^2{gbet) 
is well defined, calculable, independent of any regularisation method, and at the limit of 5 — )■ 
(or /? = Arriy^/rn^j 0) is 



M 



2(GBET) 



2e^g 



/3^0 . 



(2.19) 



By equating eq. ( 2.15| ) (in the limit /? — )• 0) and eq. (2.19) which represent the l.h.s and r.h.s of 
the GBET condition (2.17), respectively, we find 



A = -1 



(2.20) 



This value agrees with dimensional regularization 129 



30 



in the limit d — ^ 4 [see eq. (2.13)]. The 



final form of the A^2 in eq. (2.2) is 



(47r)^ my/ 



2 + 



3/3 + 3/3(2-/3)/(/3) 



(2.21) 



with /3, and /(/3) defined in eq. (2.16). 



To complete the picture there is still the coefficient M.i in eq. (2.2) to be calculated. Naive 



power counting says that this is by two powers more divergent than and, in general, un- 
determined. It can be fixed however by using quantum gauge invariance i.e., conservation of 
charge, for the C/(l)em, 



, kouM 



0, 



^2 



0, 



and thus from eq. (2.2) 



Ml = -{ki-k2)M 



2 • 



(2.22) 



(2.23) 



Eq. (2.23) is substituted to eq. (2.2) with A^2 read by eq. (2.21). This is exactly the same result 



for the VF-boson contribution to -ff — )• 77 amplitude, that has been obtained in refs. jll -14,19 
using dimensional regularisation in i?^-gauges. 

It is interesting here to note the result from the explicit algebraic manipulation oi Mi in 
the unitary gauge and check the validity of gauge invariance [eq. (2.23)]. Exactly as for 7^2, 
the condition a = b = c for the arbitrary vectors given in eq. (2.7) is crucial in reducing the 



divergence of A4i down to a logarithmic one [see expression eq. (A. 11)]. In d-dimensions the 



expression for A^i is finally independent of any arbitrary vector and, up to a proportionality 
factor, reads: 



Ml 



dx 



dy 



d'^e r 4(£ • ki)i£ ■ k2) + 2a - i)e{ki ■ k2) + (^)(4 - d) e m 



(27r)^ 



w 



+ 



{d- 



- A)3 

3 wby^ rrijj + (1 — d) x{x + y — l)m 



+ 



(£2 _ A)3 



(2.24) 



Clearly the first integral in eq. (2.24) is ill-defined in four dimensions. If however, we insist in 
doing the calculation of eq. (2.24) in d = 4 with symmetric integration, like in refs. |32l^3), 



we find that gauge invariance [eq. (2.23)] is not satisfied. This is of course unacceptable. By 
going a little bit deeper, gauge invariance is lost because of the term proportional to 4 — d in 
eq. (2.24) which vanishes when d = 4. Quite the contrary in DR, this term results in a non-zero 
contribution when mw 7^ 0, since the (log divergent) integral in front of (4 — d) contains a simple 



pole at d = 4. This changes the final result and renders eqs. (2.21), (2.24) and (2.23) consistent, 
only if A = — 1. This outcome is in agreement with ref. |19]. 

Few remarks are worth mentioning here. Had we started calculating Mi first there would 
be no possibility of defining unambiguously A without using a gauge invariant regulator: the 
Qfj^y part of the amplitude at g — )• involving Goldstone bosons [see diagrams fig. [2] is not well 



defined - an integral as the one in eq. (2.12) appears again. Another remark is that the same 
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expressions for the coefficients Aij displayed in Appendix A in the unitary gauge, appear also 
when one exploits the i?^-gauge. In the latter there are in addition ^-dependent terms fl9] that 
vanish in the end from unphysical scalar contributions. Therefore, the logarithmic ambiguity in 
eq. (2.12), found here in the unitary gauge, is similar in every other gauge. 



3 Four Dimensional Regularization (FDR) 



So far we have proposed a regularization scheme which is four-dimensional and uses the basic 
symmetries and underlying physics of the SM. However, in more complicated models or observ- 
ables with more parameters to adjust, such a scheme can become cumbersome. For example, it 
is not always obvious which physics argument will fix undefined integrals. 

Very recently, R. Pittau (31] proposed a scheme that is fairly easy to handle and, to the best of 
our knoweledge, is the closest to four dimensional calculations, thereby coined four-dimensional 
regularisation/renormalization scheme or just FDR. According to this scheme, infinite bubble 
graph contributions, i.e., large loop momenta contributions that do not depend upon external 
momenta, are absorbed into the shift of the vacuum while the remaining finite corrections are 
calculable in four-dimensions in addition to being Lorentz and gauge invariant. 

We have applied FDR into the calculation of the — t- 77 amplitude and found agreement 
with our physics approach and with DR results. In FDR one introduces an arbitrary scale /i 
which is considered to be much smaller than internal momenta and particle masses in loops. Self 
contracted loop momenta quantities like become P = i"^ — iJ,'^, while for gauge invariance to 
hold, vector momenta, p^, remain untouched. For example the integral of eq. (2.12) becomes. 



id' 



(3.25) 



where D = [P — A) and [d'^t] stands for integration over dropping all divergent terms from 
the integrand (see below) and taking the limit ;U — )• 0. In going from l.h.s to r.h.s of eq. (3.25 ) the 
symmetry property i^li, = g^yP/^ has been used in four dimensions. Then, using the partial 
fractions identity. 



1 



+ A 



1 



+ 



1 



1 



(3.26) 



the term in square bracket is recognised as divergent and therefore removed, and integrating the 
r.h.s of eq. (3.25) over {d'^P\ one obtains 



[d' 



-A lim 



1 



+ 



1 



+ 



1 



(3.27) 



i.e., exactly the same result as in DR which eventually leads to A = — 1 consistent with gauge 
invariance and GBET. What in fact FDR scheme does, is to restate the correct DR answer 
through the regulator /i2 keeping eq. (2.12) correct in d = 4. We therefore understand that the 



constant (/3-independent) term of eq. (2.21 ) in FDR arises from the fact that the arbitrary scale, 
/i^, must disappear from physical observables. 
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4 Discussion 



It is evident that our calculation for the amplitude incorporates two physical inputs: one is the 
conservation of charge and the other is the equivalence theorem. They are both direct conse- 
quences of the gauge invariance of the underlying physical theory. The first one is experimentally 
indisputable while the second one is theoreticaj^ and has been proven in ref. 42 that is valid 



in any spontaneously broken renormalizable theory, like for example the SM. One may think 
however that there is a loophole in our use of this second argument: so far, and, to our knowl- 
edge, the replacement of the VF-bosons with Goldstone bosons at high energies has been proven 



to be valid only for external I^-bosons |38 , 43 , 44 and not for internal ones which is the case 



exploited here. Although it has been tested in several phenomenological examples [45] , a formal, 
to all orders, proof is still missing. Although this may be true, it is difficult to argue against the 
validity of decoupling limit 5 — )• (with fixed v.e.v and Higgs quartic coupling) discussed in the 



paragraph above eq. (2.18). 



Is there another physics context from which one can define A? One possibility is to exploit the 



low energy Higgs theorem 11,46-48 instead. Although this may serve as a consistency check, 
and indeed is compatible with A = —1, we cannot use it to define A. The reason here is threefold: 
first, when treating the Higgs field as an external background field with zero momentum one 
needs to take partial derivative w.r.t mw of the 2-point photon vacuum polarization amplitude, 
n^^(g^). The later, is notoriously difficult, if meaningful, to be calculated in the unitary gauge. 



Second, according to ref. 11 , we know that to the lowest order in weak coupling, the amplitude 
for the process {'yj\H) is proportional to (77|O/i'|0) where @li = 2ml^W^W~+... is the improved 
energy momentum tensor ^49j. However, the calculation of (7710/^*10) goes through the same 
steps as for the calculation for the — >• 77 amplitude and therefore involves the same ambiguity 
for calculating A. Third, one could examine the H^-contribution to i/ — t- 77 within the dispersion 
relation approach. It can be shown [SO] that the non-vanishing limit at gw — )- is due to a finite 
subtraction induced by the corresponding trace anomaly j51j. However, in order to calculate 
unambigiously this finite piece, one has to make full use of a (physical) boundary condition 
of the theory. In conclusion, an unambiguous calculation of the — )■ 77 amplitude must be 
performed within a physical setup of gauge invariance and the GBET. 

As a final remark, suppose that we did not know DR and wanted to calculate a certain 
observable in 4-dimensions. In this observable we encounter singularities i.e., undefined and 
undetermined integrals. Then we use physics arguments to fix these ambiguities. However, we 
can always question whether we are using the right physics set up or not. In that sense the final 
judgement should come from the experiment. Therefore it may be not academic to ask whether 
LHC could see the difference between A = — 1 and A = 0? Setting the SM Higgs mass rriH = 125 
GeV, and including the top-loop contribution, we find 

Bi{H ^ 77, A = 0) , , 

— ^- '- « 0.46 . 4.28) 

Br(iJ ^ 77, A = -1) ^ ' 

This is certainly within LHC's sensitivity for 14 TeV cm energy and luminosity of 30 (see 



This is not entirely correct. There is of course the high energy behaviour of e e — >■ W W found at LEP [41] 
consistent with the GBET. 
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for example Fig. 3 in ref. 152 



„ valiiP Br(H^77, (exp)) 
^ ^^^^^ Br(//^77,A=-1) 



1.6 ±0.3 



53 



In fact, the recent observation by LHC experiments [6|7] indicates 
which highly disfavours the case A 



by almost four 

standard deviations. We can turn this around and state that this is an indirect hint towards the 
validity of the equivalence theorem. 



5 Conclusions 



In this work we review the W-g&uge boson loop contribution to the H ^ jj amplitude in the 
unitary gauge. Our objective is to fix intermediate step indeterminacies arising from divergent 
diagrams by making full use of physics at d = 4 much in the same way as in the calculation of 
the chiral anomaly triangle. 

We anticipate a finite result for the loop induced H — )• 77-amplitude in the renormalizable 
SM. Therefore the amplitude has to be independent of any shifting momentum variables we 
have originally introduced. But finite or even log divergent integrals are independent of these 
vectors, so the vectors have to be accompanied only by infinite contributions, if at all. Therefore, 
infinities and arbitrary vectors are eliminated altogether by a certain combination among them 



[see eq. (2.7)]. 



The whole calculation in the unitary gauge boils down to a logarithmically divergent integral 
(2.10). We find that, this integral results in two different values depending on whether d — 4 
or d = 4. This is due to a surface term remaining at the exact d = 4 case after the part-by- 



part integration in d-dimensions [see Appendix B . To proceed, we identify this integral with an 



undefined parameter A [see eq. (2.12)]. This parameter is then fixed unambiguously by assuming 



the validity of the Goldstone Boson Equivalence Theorem (GBET). Its value is consistent with 
DR in the limit d — > 4. 

In our calculation we are very careful not to perform shifting of integration variables for highly 
divergent integrals by introducing three arbitrary momentum variable shift vectors straight from 
the beginning. Divergencies and arbitrariness from these unknown vectors are altogether re- 



moved, leaving behind a log-like divergent integral in A^2 of eq. (2.2). This is defined by a 



physical input taken from the GBET and is connected to TWi by electromagnetic charge conser- 
vation. 

As noted many times in the text, the key point towards deriving the unambiguous amplitude 
for — )• 77 in the unitary gauge is the limit of vanishing gauge couplings; this is an aspect of 
GBET [eq. (2.17)]. In this limit, the Goldstone boson loop contribution to the coefficient A^2 is 



finite, independent of any regularisation scheme. 

We also saw that DR (FDR) , a regularisation scheme introduced to maintain Ward Identities 
at intermediate steps of a calculation, supports the GBET in the limit d — >• 4 (d = 4). On the 
contrary, we find that, performing the integrals in d = 4 with symmetric integration is not a 



good choice because it leads to the violation of gauge invariance [see eq. (2.23) and the discussion 
below] . The main reason is due to surface terms that are developed in exactly d = 4 dimensions 



[see discussion below eq. (2.10) and Appendix B| . The latter are axiomatically discarded in 
DR 129.(30]. Another reason is the appearance of the (d — 4)-term in the numerator of eq. (2.24). 
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In conclusion, the four-dimensional calculation of — )• 77 amplitude in the unitary gauge 
is ambiguous without introduction of a physics input beyond gauge invariance. As we have 
demonstrated, this physics, which uniquely defines the amplitude, arises from the Goldstone 
Boson Equivalence Theorem (GBET). This effectively proves that GBET comprises an additional 
important pillar of the Standard Model dynamics. 

Acknowledgments: We are grateful to Adrian Signer for illuminating discussions and crit- 
ical comments on the manuscript. We would also like to thank P. Kanti, S. Martin and K. Tam- 
vakis for discussions and comments. We also thank C. Coriano for bringing to our attention 
ref. fSO]. This research Project is co-financed by the European Union - European Social Fund 
(ESF) and National Sources, in the framework of the program "ARISTEIA" of the "Opera- 
tional Program Education and Lifelong Learning" of the National Strategic Reference Frame- 
work (NSRF) 2007-2013. K.S. acknowledges full financial support from Greek State Scholarships 
Foundation (I.K.Y). 



Appendix A 



We append here the integrand expressions for the coefficients Aij in eq. (2.3). The corresponding 
formula for is quite long and is not included here. It can be provided by the authors upon 
request. Note that the number of dimensions d has been kept arbitrary throughout and on-shell 
conditions for the external particles have been imposed. 

1 



A' 



21 



+ a)2 - m'^][{p + a - ki^ - m'L][{p + a - ki - ^2)^ 



m 



Wi 



(4 d - 6) + 



+ 



+ - 



m 



W 



3 {p + a) ■ {p + a) - 5 {p + a) ■ ki - {p + a) ■ k2 + 2 m^j 
Hp + a)-{p + a)f + 3 ((p + a) • A;i)((p + a)-{p + a))-2 {{p + a) ■ kif - 
-2 {{p + a) ■ ki){{p + a) ■ k2) + {{p + a)-{p + a)){{p + a) ■ k2 



A22 



[{p + 6)2 - m|j^] [{p + b- k2)^ - rr^y] [{p + h-ki- ^2)^ - m 



(4 d - 6) + 



+ - 



m 



w 



3 (p + b) ■ (p + b) -5 (p + b) ■ k2- (p + b) ■ ki + 2 mjj 

Hp + b)-{p + b)f + 3 ((p + 6) • k2){{p + b)-{p + b))-2 Hp + b) ■ k2f 
-2 Hp + b) ■ ki)Hp + a) ■ k2) + Hp + h)-{p + b))Hp + b) ■ k^) 



+ 



(A.l) 



(A.2) 
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A23 



[{p + c)2 - + c-ki- k2f - m^] 



4 + 



m 



w 



{{p + c)-{p + c)) + {p + c)-ki + {p + c)-k2 



A31 



[{p + a)2 - m'^][{p + a - A;i)2 - m'^][{p + 0-/^-/^2) 



m 



(7 - 4 d) - 



4 (p + a) • (p + a) - 7 (p + a) • fci + 3 + a) • A;2 + 
((p + a) • (p + a))2 - 3 ((p + a) • A;i)((p + o) • (p + a)) + 2 ((p + a) • fci)^ + 
+2 ((p + a) • ki){{p + a) ■ fca) - ((p + a) • (p + a))((p + a) • /c2 



A32 = 



[{p + b)^-ml,][{p + b-k2r- ml.] [{p + b-ki- k2Y - ml,] 
rriw + - (p + b) ■ (p + b) + 6 (p + b) ■ k2 



A: 



•33 



2 - 



m 



w i 



[{p + c)2 - m^^] [(p + c - /ci - ^2)2 - m^] 

(p + C) • (p + C) - (p + C) • fcl - (p + C) • ^2 



A4I = 



[(p + a)2 - m^] [(p + a - A;i)2 - m^] [(p + a - fei - ^2)^ - m^] 
"i^w + — (i? + a) • (p + a) + 6 (p + a) • fei 
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+- 



m 



w 



"^^^ [{P + hf - m^y] [{p + b-k2)^ - m"^] [{p + b - ki - k2)^ - m^^] 
(7 - 4 d) - A (p + b) ■ (p + b) -7 (p + b) ■ k2 + 3 {p + b) ■ ki 

{{p + b) ■ {p + b)f -3 {{p + b) ■ k2){{p + b) ■ (p + b)) + 2 {{p + b) ■ k2f + 

+2 {{p + b) ■ ki){{p + b) ■ k2) - Hp + b)-{p + bMp + b) ■ ki) 

1 



(A.8) 
(A.9) 



51 



[{p + a)2 - m^] [{p + a- kiY - ni^] [{p + a - ki - k2Y - m^] 



5m^ + 



+ 



3 (p + a) -{p + a) -2 (p + a)- ki 
1 



+ 



[{p + 6)2 - mh] [{p + b- k2)'^ - m'L] [{p + b - ki - k2Y - m 



5 myi/ + 



3 (p + 6) • (p + 6) - 2 (p + 6) • A:2 
2 

\[p + c)2 - m^] [(p + c - A;i - /C2)^ - m^] 



(A.IO) 



It is straightforward, but long and tedious, to show that after implementing the condition (2.7) 
to coefficients in eqs.(A.l)-(A.10) we arrive at eq. (2.8) which is at the most logarithmically 
divergent. 



For complementarity reasons, it is useful in deriving eq. (2.24) to present the expression for 
ien 

^11 



the coefficient after the imposition of the arbitrary vector relation eq. (2.7). 

1 



[{p + a)2 - m^] [{p + a-kiY - m^] [{p + a - ki - k2)' 



m 



{p + a- kif - ] {I - d) rn^ + 



+ A[{p + a) ■ ki] [{p + a) ■ /ca] - [3m^ + {p + ay] mij ) + 



[{p + a)2 - rn^] [{p + a - k2Y - m1^] [{p + a - ki - k2Y 



m 



{p + a- k2f' - rn^j (1 - d) + 
+ A[{p + a)- ki] [{p + a) ■ /C2] - [3m^ + {p + af\ m% 



(A.ll) 



This integrand expression, under J d'^p, is obviously at the most logarithmically divergent. 
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Appendix B Dimensional Regular izat ion and the surface term 



We would like to examine the surface terms arising in d = 4 when calculating the integral on 



the l.h.s of eq. (2.12). This integral after Wick rotation into Euclidean space, reads 



Qui [I ^ dflU 4 ly 



(B.l) 



where i = Ie and drop for clarity the subscript E from now on. We follow very closely 't Hooft 
and Veltman's seminal paper in ref. |29| . In our calculation for a physical process we should 
notice first that i^,£u are strictly 4- vectors since they are contracted with physical external 
momenta or /cj'g- the other hand, the loop momentum £ in has components in all, 
d = 2uj, dimensions. We write £ as a sum of a vector which has non-zero components in 
dimensions 0, 1, 2, 3 and a vector £± which has nonzero components in (2a; — 4)-dimensions, 



£ = £\i+ £± 



With this definition, the integral (B.l) reduces to 



i / d^'^£ 



'2lu p ^1 9lJ.v 



(^2 + A)3 ' 



(B.2) 



(B.3) 



where the £11 components in the numerator of (B.l) vanish thanks to symmetric integration 



formula, — ?■ \£'^\^g^'^ ■ In order not to carry the g^^y in all formulae below we just concentrate 
on the integral 



+ 4 + A)3 



Integrating over the extra dimensional solid angle d02aj-4 we arrive at 



r(w - 2) 



dH 



dL 



L 



(£2 + L2 + A)s 



(B.4) 



(B.5) 



where r(x) is the Euler F-function and L is the length of the £± vector. This integral is UV 
divergent for uj >2 and IR divergent for a; < 1. Therefore, the region of convergence, 1 < w < 2, 
is finite but it does not yet include the point uj = 2. In order to enlarge the region of convergence 
to include u = 2 one has to analytically continue I by inserting the identity. 



1 



d£ 



II M 



d£ 



+ 



dL 

dL 



(B.6) 



in (B.5). After integrating by parts in the region of convergence, rewriting the r.h.s in terms of 



I from eq. (B.5) and keeping only, potentially, non- vanishing surface terms, we arrive at 

-2A 



T(4 



UJ 



^11 M 



6 in'' 



Tiu - 1) 



dL 



L 



2lu-3 



{£l + L2 + A)4 



(B.7) 
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where the first integral is over the Euclidean spatial components of a 4- vector on a three-sphere. 
The surface integral converges in 1 < cj < 2 while the other in 1 < a; < 3. By taking the surface 
integral on a three-sphere with radius R and eventually taking the limit i? — ?■ oo we find 



/ 



27r^ hm R 



.2^-4 



(B.8) 



which now converges in the region u < 2, that is, it includes the point lo = 2. For lo < 2 this 
surface term vanishes while for ui = 2 there is a finite piece, 27r^, remaining. This is exactly the 
term that spoils gauge invariance and the equivalence theorem. In DR this term is axiomatically 
absent - the shifting of integral momenta is among DR's main properties. 



Turning into the second integral of eq. (B.7) we note first that the region of convergence 
includes now oj = 2. It gives, 



-2a;- 3 



dL 



tt' Tjoo 
12 



l)r(3 



(B.9) 



By placing eqs. (B.8) and (B.9) into eq. (B.7) we finally arrive at eq. (2.12). 
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